An Equation of State of a Carbon-Fibre Epoxy Composite under Shock Loading 
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An anisotropic equation of state (EOS) is proposed for the accurate extrapolation of high-pressure 
shock Hugoniot (anisotropic and isotropic) states to other thermodynamic (anisotropic and isotropic) 
states for a shocked carbon-fibre epoxy composite (CFC) of any symmetry. The proposed EOS, 
using a generalised decomposition of a stress tensor [Int. J. Plasticity 24, 140 (2008)], represents 
a mathematical and physical generalisation of the Mie-Griineisen EOS for isotropic material and 
reduces to this equation in the limit of isotropy. Although a linear relation between the generalised 
anisotropic bulk shock velocity Uf and particle velocity u v was adequate in the through-thickness 
orientation, damage softening process produces discontinuities both in value and slope in the Vf-u p 
relation. Therefore, the two-wave structure (non-linear anisotropic and isotropic elastic waves) that 
accompanies damage softening process was proposed for describing CFC behaviour under shock 
loading. The linear relationship U s -u p over the range of measurements corresponding to non-linear 
anisotropic elastic wave shows a value of Co (the intercept of the V^-u p curve) that is in the range 
between first and second generalised anisotropic bulk speed of sound [Eur. Phys. J. B 64, 159 
(2008)]. An analytical calculation showed that Hugoniot Stress Levels (HELs) in different directions 
for a CFC composite subject to the two-wave structure (non-linear anisotropic elastic and isotropic 
elastic waves) agree with experimental measurements at low and at high shock intensities. The 
results are presented, discussed and future studies are outlined. 
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INTRODUCTION 

Investigation of anisotropic composite materials (e.g., 
CFC materials) behaviour has found significant interest 
in the research community due to the widespread ap- 
plication of anisotropic composite materials in aerospace 
and civil engineering problems. Modern high-resolution 
methods for monitoring the stress and particle velocity 
histories in shock waves and equipment have been cre- 
ated [J 0, d, 0> S] • A common technique for the study 
of material behaviour under shock loading is the planar 
plate impact test (one-dimensional shock wave propaga- 
tion) . Shock wave experiments have frequently provided 
the motivation for the construction of material consti- 
tutive relations and have been the principal means for 
determiningmaterial parameters for some of these rela- 
tions [1, 0, H, & Ell- For example, Dandekar et al. [TlT ] 
investigated the equation of state of a glass fibre epoxy 
composite, in terms of the shock stress, shock velocity 
U s (through-thickness orientation) and particle velocity 
Up (i.e. the velocity of material flow behind the shock 
front). Their results indicated that there was a linear 
relationship between shock and particle velocity. This 
type of behaviour is typical of a wide range of materi- 
als, including metals H Q and some pol yme rs, includ- 
ing epoxy resins 12|, |l3[ , and composite [14j , including 
carbon fibre epoxy composite [l5| and glass fibre epoxy 
composite 16]. A linear U s -u p relationship shows that in 
the through thickness orientation, this class of compos- 
ite displays fairly typical experimental data. However, 
in spite of a perfectly adequate general understanding, 
experimental methodology, and theory, material models 
do not agree in detail, especially for anisotropic compos- 



ite materials. For many years, it has been assumed that 
the response of c omp osite materials to shock loading is 
[lH fig ], and only recently has anisotropy 
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isotropic 

in the shock response of anisotropic materials (e.g., com- 
posite materials ) a ttracted the attention of researchers 

HEim, QUI III- 

In this paper, a macroscopic continuum framework is 
chosen to modify the existing methodology for accurate 
extrapolation of high-pressure shock Hugoniot states to 
other thermodynamic states for a shocked isotropic con- 
tinuum. The composite materials response under shock 
loading leads to a nonlinear behavior (i.e. large compres- 
sions), therefore, an equation of state (EOS) is required 



[24 I2BI |26| ]. To address this issue, an anisotropic 
equation of state similar to 24, 25, [5^] is proposed for 



accurate extrapolation of high-pressure shock Hugoniot 
(anisotropic and isotropic) states to other thermody- 
namic (anisotropic and isotropic) states for a shocked 
carbon-fibre epoxy composite of any symmetry. 



AN ANISOTROPIC EQUATION OF STATE 



The proposed equation of state, usin g[ th e generalised 



decomposition of a stress tensor [3, [25, 2f| 27, 28 1, rep- 
resents a mathematical and physical generalisation of the 
Mie-Griineisen equation of state for isotropic material 
and reduces to this equation in the limit of isotropy. The 
generalised decomposition of the stress tensor is sum- 
marised below. The generalised decomposition frame- 
work will provide a useful point of construction of an 
anisotropic equation of state. 
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Generalised decomposition of the stress tensor: a-(3 
decomposition 

The definition of pressure in the case of an anisotropic 
solids should be the result of stating that the "pressure" 
term should only produce a change of scale, i.e. isotropic 
state of strain O, HI, [H, [2]], Hj| . The generalised de- 



composition of the stress tensor aij is defined as: 



P : S = 



OLijSij = 0, o-jj = -p*aij + §ij, (1) 



where P = p*a.ij is the generalised spherical part of the 
stress tensor, S = Sij is the generalised deviatoric stress 
tensor, p* is the total generalised "pressure" and ctij is 
the first generalisation of the Kronecker delta symbol. 
The constructive definition of the tensor cnj is based 
on the fact that the stress tensor paij is induced in the 
anisotropic medium by the applied isotropic strain tensor 
yd i:) , i.e. 



HjklOkl-, 



P = -K c e v , 



(2) 



where p is the pressure, s v is the volumetric strain, Sij 
is the Kronecker delta symbol (unit tensor), Cijki is the 
elastic stiffness matrix and Kc is the first generalised 
bulk modulus. The expressions for the a,ij components 
and Kc are presented below. Throughout, the contrac- 
tion by repeating indexes is assumed. Using (Q}, the fol- 
lowing expression for total generalised "pressure" p* can 
be obtained: 



ctkiaki 



1 



2 a ij ' 



(3) 



where K$ represents the second generalised bulk mod- 
ulus. The solution of equations (J6]) in terms of the j3ij 
components and an expression for K$ are presented be- 
low. Equations (|5|) define the correct generalised "pres- 
sure" for the elastic regime. Note that the generalised 
decomposition of the stress tensor can be applied for 
all anisotropic solids of any symmetry and represents a 
mathematically consistent generalisation of the conven- 
tional isotropic case. The procedure of construction for 
the tensor oiki has been defined in [3, [3, Hi" 
The elements of the tensor otki are 



27, 2b_ 
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(8) 

where dj is the elastic stiffness matrix (written in Voigt 
notation) . The elements of the tensor j3ki are 

fti = (E Jk ^j ■ 3Ks > & 2 = (E Jfc2 ) ' 3 ^ 5 ' 

I3 33 = (E Jfc3 J - 3Ks ' fa Pa = 3 ' 



(9) 



vfc=i 



where ||a|| = ay-a,j = a\ Y + a^ 2 + "33- Finally, the 
expression for the generalised deviatoric part of the stress 
tensor can be rewritten in the following form: 



Si 



1 



r^kiaki- 



For anisotropic materials, the totalgeneralised ' 
sure" p* has been expressed 24, 25, 26, 27, 28| as: 



pres- 



P + P > P 



Pij a ij s 

a > " 
PijCtij 



Pij Sjj 
PklCtkl 



(5) 



where p is the pressure related to the volumetric defor- 
mation ([2]), p s is the pressure related to the generalised 
deviatoric stress and Pij is the second generalisation of 
the Kronecker delta symbol. The constructive definition 
of the tensor /5y is based on the fact that stress tensor 
p/3ij is independent of the stress tensor Cijki&kU i- e - their 
contraction product is zero for any deviatoric strain ten- 
sor e^i, where Gyjy is the elastic stiffness matrix. The 
following relation describes the functional definition of 
the second order material tensor P^: 



Ks 



\ 




(4) where Jij are elements of the compliance matrix (writ- 
ten in Voigt notation). In the limit of isotropy, the pro- 
posed generalisation returns to the traditional classical 
case where tensors otij and fyj equal Sij and equations 
([5]) take the form: 



c x = -n a kk 
OklOkl O 
1 



PklOkl 
Pij Sij 
Pkl^kl 
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(11) 



0. 



Cfcfc 



where p = — — is the classical hydrostatic pressure. 

Also, parameters Kc and K$ were considered as the first 
and the second generalised bulk moduli. In the limit 
of isotropy, they reduce to the well-know expression for 
conventional bulk modulus. 



PijCijki = SKsSki, 



(6) 
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or 



Pi Pi 



FIG. 1. Schematic representation of a-decomposition of the 
stress tensor (Haigh-Westergaard stress space). 
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, . 70 + ap p 
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1 + M Po 



(14) 
(15) 



where 70 is the initial Griineisen gamma, a is the first 
order volume correction to 70, Pg is the anisotropic gen- 
eralised bulk Hugoniot pressure, E is the internal energy 

per initial density [ E = — ) , a is the relative change of 
V PoJ 

volume, v is the specific volume and Sf, S^, S3 are the 
intercept of the cubic U^-u p curve 0]: 



The geometrical representation of generalized decom- 
position of the stress can be shown in the principal 
stress space (Haigh-Westergaard stress space) for a- 
decomposition of the stress tensor and /^-decomposition 
of the stress tensor in the principal strain space. FigureQ] 
shows schematic representation of a-decomposition of the 
stress tensor, where a-direction is described by the tensor 
ctij and (5-direction is described by the Kronecker's delta 
tensor 5%j. Therefore, p$ describes the hydrostatic stress 
(isotropic stress), p* ^ p§ describes the total generalized 
hydrostatic stress (or anisotropic total generalized hydro- 
static stress) , and p 7^ p* ^ pg is the generalized pressure 
related to an equation of state (EOS). The angle between 
a-direction and ^-direction is described by the variable 

-0 which can be obtained from cos \j) = — — ^ — . 

Similar representation can be shown for /3-decomposition 
of the stress tensor (Haigh-Westergaard stress space) . 



An anisotropic EOS 

The equations j5]) define the correct generalised "pres- 
sure" for the elastic regime: 



P 



P- 



~P OLij 



St. 



(12) 



Further, to provide an appropriate description for gen- 
eral anisotropic materials behavior at high pressures, the 
pressure p related to the volumetric deformations is de- 
scribed by the Mie- Griineisen equation of state pressure 



V 



EOS. 



P = P 



EOS 



f(p,e)=Pj 



H 



pTe (13) 
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where Uf is the generalised anisotropic bulk shock ve- 
locity and Cq is the velocity curve intercept. The pro- 
posed generalised "pressure" also correctly describes the 
medium behavior at small volumetric strains. To be con- 
sistent with the definition of the isotropic bulk speed of 
sound, the following definitions of the first cj and the 
second en bulk speed of sound for anisotropic medium 
are assumed 1251: 



ci 




en = 




(17) 



where the generalised bulk moduli Kc, Ks are defined 
according to (jHJl and (fl"U)) respectively. Parameters Cq e 
[cii,cj], S A , S%, S3, 70, a represent material properties 
which define its EOS (fill). 



THE BEHAVIOUR OF A CARBON-FIBRE 
EPOXY COMPOSITE UNDER SHOCK LOADING 

The purpose of this section is to display the ac- 
curate extrapolation of experimental Hugoniot states 
[23| behind the shock wave (thermodynamic states) 
to high-pressure shock Hugoniot Stress Levels (HSLs) 
for a shocked carbon-fibre epoxy composite (initially 
anisotropic CFC) using the anisotropic EOS proposed 
above. 



Description of Experiment 

The work discussed below concerns the shock response 
of a carbon-fibre epoxy composite. This is done by 
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the technique of plate impact, whereby a flat plate of 
constant thickness and of known material (for instance 
aluminium alloy, or copper) is impacted onto a target 
plate made from the test material. The flyer plates are 
launched using a 50 mm bore, 5 m long single stage 
gas gun. The plate impact test was done at the De- 
fence Academy of the United Kingdom by Millett et al. 
[23| using samples of a carbon-fibre composite (CFC) of 
thicknesses 3.8 mm. The impact axises were (a) nor- 
mal to the plane of the fibres and (b) parallel to the 
plane of the fibres. On impact, a planar shock front 
starts propagating into the target. The shock propaga- 
tion in the target is monitored using manganin stress 
gauges, placed at different locations within the target as- 
sembly. A schematic of the target assembly and gauge 
placement is shown in Figure [2] A manganin stress gauge 
was supported on the back of the specimen plate with a 
12 mm block of polymethylmethacrylate (PMMA). Also, 
the gauge was backed into the PMMA by approximately 
1.5 mm PMMA offset block to act as extra protection for 
the gauge. A second gauge (the mm position) was sup- 
ported on the front of the target assembly with a 1 mm 
plate of aluminium alloy 6082-T6. Shock stresses were 
induced with dural flyer plates impacting with the dif- 
ferent velocities, using a single stage gas gun (23| . The 
results from the stress gauges were converted to material 
(Target) values <jm, using the shock impedances of the 
target At and PMMA Ap, via the well-known relation 



A T + A P 
° M = -^Ap— ap > 



(18) 



where ap is the stress gauges values. Furthermore, tem- 
porary spacing (St) between the tracers in combination 
with the specimen thickness (Sw) were used to obtain 
the shock velocities in the longitudinal (through thick- 
ness orientation [7~ = Sw/5t and along the fibre 0° ori- 
entation Uf = 6w/8t, taking into account both the offset 



FIG. 3. Microstructure of carbon-fibre epoxy composite. 



distance of the gauge within the PMMA and the known 
shock response of PMMA d, [13|. 



Materials 

The fibres in the carbon-fibre epoxy composite were 
Hexcel 5HS in a woven lay up of orientation [0/90, ±45]4. 
The arial weight was 370 g/m 2 . The resin was an epoxy, 
Hexcel RTM6. The composite was cured at 180° C un- 
der a pressure of 100 psi (6.7 atm) for 1 hour 40 minutes. 
The microstructure is shown in Fig. ([3]). The longitudi- 
nal sound speed in the through-thickness orientation Cl 
was 3020 m/s, and the ambient density p was 1.50 g/cm 3 . 
The equivalent material properties of the CFC compos- 
ite plate were chosen to match the layer macromechan- 
ical properties for the layup [0/90, ±45]4 only and the 
longitudinal sound speed in the through-thickness orien- 
tation. A description of their numerical values for the 
selected Carbon Fibre Composite (CFC) is shown in Ta- 
ble HI The acoustic longitudinal speed of sound in the 
through-thickness orientation is: 



Cr 



'IBol 

.P d£ Lj p=po 



(19) 



where <7l is the longitudinal stress, El is the longitu- 
dinal strain. Using the generalised Hook law for or- 
thotropic materials and an uniaxial acoustic strain state 
(through thickness orientation), i.e. El = £33 = £zz 
(£ij ^ i,j = 3 only), the longitudinal speed of sound 
(through thickness orientation) fl9]) takes the form: 



Cr 



(1 - 1/121/21) E 3 



(20) 
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TABLE I. Material properties of undamaged (CFC) Carbon 
Fibre Composite (11-direction corresponds to the fill direc- 
tion, 22-direction to the warp direction and 33-direction the 
through-thickness orientation). 



Parameter 


Description 


CFC 


p$ [kg/m 3 ' 


Density 


1500.0 


E 1 [GPa] 


Elastic modulus (22 direction) 


68.467 


E 2 [GPa] 


Elastic modulus (11 direction) 


66.537 


E s [GPa] 


Elastic modulus (33 direction) 


13.678 


K c [GPa] 


Firtst generalised bulk modulus 


19.436 


K s [GPa] 


Second generalised bulk modulus 7.6902 


C L [m/s] 


Sound speed (11 direction) 


6762.0 


Cl [m/s] 


Sound speed (22 direction) 


6666.0 


Cl [m/s] 


Sound speed (33 direction) 


3020.0 


V%\ 


Poisson's ratio 


0.0400 


VM 


Poisson's ratio 


0.0045 


VZ2 


Poisson's ratio 


0.0044 


„A 

a n 


Tensor a (22 direction) 


1.2290 


&22 


Tensor a (11 direction) 


1.1956 


„ A 
«33 


Tensor a (33 direction) 


0.2454 


A A i 


Tensor j3 (22 direction) 


0.3155 


022 


Tensor j3 (11 direction) 


0.3254 


$53 


Tensor (3 (33 direction) 


1.6717 



where A = (1 - f 121*21 - ^23^32 - ^31^13 - 2^21^32^13); 
Ei, E2, E3 are Young's moduli and are Poisson ra- 
tios. The measured longitudinal speed of sound was 
Cl = 3020 ±30 m/s for the initial material density 
p = 1500 ± 10 kg/m 3 Using flU and elastic mate- 
rial properties in Table [H the longitudinal speed of sound 
is calculated to be Cl — 3020 m/s, which is in exact 
agreement with the measured value of 3020 ± 30 m/s. 
Note that we have a six elastic unknown material prop- 
erties listed in Table [l] which directly affect shock waves 
propagation under uniaxial strain state. However, only 
one longitudinal sound speed in the through-thickness 
orientation, depending on four unknown elastic proper- 
ties, was experimentally measured (23[. Hence, the equiv- 
alent material properties given in Table U have not been 
finalised yet for the selected CFC material. 



Results and discussion 

Initially, the equations of state, in terms of the shock 
and particle velocities, arc examined. Although a linear 



relationship in the U, 



A.L.F 



i p plane is adequate for sim- 



ple materials and some anisotropic metals [24, 25, 26, 29], 
phase changes under shock loading (e.g., damage soften- 
ing) produce discontinuities both in value and slope in 
the U^' L ' F -u p relation, where Uf' L ' F (i.e., Uf, C/ S L , Uf) 
represent the generalised shock velocities in the directions 



of anisotropic bulk, longitudinal through-thickness, and 
along the fibre 0° orientation, respectively. These dis- 
continuities are usually caused by the two- wave structure 
that accompanies most phase changes, as well as damage 
softening. It was pointed out by Bethe [3f| that, for sta- 
ble shock waves, the shock velocity must increase with 
pressure. This means that if the shock velocity should 
decrease with pressure, then the shock front would break 
up into two or more waves, or possibly one wave with a 
continuously smeared front. 

Note that the experimental data for the through thick- 
ness orientation (see, Fig. Q) can be fitted by a linear 
relation, and there is no explicit evidence of the shock 
front breaking up; however, the analysis of the exper- 
imental data for selected CFC material [23[ (see, Fig. 
(HI) shows that the shock velocity along the fibre 0° 
orientation decreases with pressure - therefore, a two- 
wave structure is proposed for describing the experimen- 
tal data. Additionally, a comparison of the equations of 
state in terms of Hugoniot Stress Levels (HSLs) for the 
single wave and two-wave structures is performed, in or- 
der to demonstrate the accuracy of the two-wave struc- 
ture methodology. As a result, the experimental data 
[23I ] for longitudinal (through thickness) orientation and 
along the fibre 0° orientation shock wave propagation in 
the selected CFC material has been fitted by straight 
lines (for single wave and two- wave structures), that is, 
using relationships for U^-u p and U F -u p that are linear 
(see, Figs. Q, ©)• 
For a single wave structure: 



a) =c^ + S 



l Up, , 



st > 0, 

if 



V u 



b) U s = c + Si Up, Oq, S( > 0, V Up, 
and, for a two-wave structure: 



(21) 



d) ut 



AF 



-Sf L 
S[ L u 

jAF„ 



Up, Oq 



Pi °0 1 



. Sf L > 0, 
S[ L > 0, 

AF 



v 'S Ui , 



> u 



JF 



St"up,4 F ,St" >0,u p <u}, 



S[ F u p , Cq F , S F > 0, 



Up > u* f , 



, (22) 

where u\ , u*j are the transition particle velocities oriented 
through the thickness and along the fibre 0°, respectively. 
Note that, as the severity of the shock increases, the 
Hugoniot Stress Levels (HSLs) of the two orientations 
converge. This fact demonstrates that the selected CFC 
material shows isotropic behaviour at high shock inten- 
sities, and can be described as an isotropic mixture of 
epoxy bunder and fractured fibres. Hence, iij and u*j can 
be treated as transition particle velocities from the struc- 
tured anisotropic material to the unstructured isotropic 
material. From the experimental data 23j, the following 
data is defined: uf = 179.5 m/s and ity = 333.0 m/s 



'7 



00, a 



(see Table HTT|) . Note that, in the limit as 
two-wave structure methodology reduces to a single wave 
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structure methodology. This property will be used to ob- 
tain EOS data for a single wave structure using the equa- 
tions that apply to a two-wave structure. Therefore, the 
interpolation algorithm and equations below have been 
written for the two-wave structure only. Also, there is a 
degree of scatter within the experimental data; however, 
no obvious variation between shock velocities Ug, Uf 
and specimen thickness has been observed [Hj]. Tradi- 
tionally, the longitudinal (through thickness) orientation 
is tabulated better. Therefore, this direction was chosen 
for an accurate extrapolation of experimental HSLs to 
EOS data. 

Using a relation XJ^-u v and the Rankine-Hugoniot 
equation expressing the conservation law of mass for a 
coordinate system in which the material in front of the 



shock wave is at rest, the compressibility factor ( — 
calculated using, 



p{U L s 



U 



(23) 



Shock wave loading deals with large finite strains. Hence, 
the Hencky strain (through thickness orientation), e H , 
can be used to describe a uniaxial strain state for a com- 
pression under shock loading [3lj |: 



fc 33 


-(f)- 


.-(! 


e H 


— fc 33' £ 11 — 


F H - 
fc 22 ^ 


"33 


- 1 F H d H 

— ^ > "11 


- d H 

— "22 



u 



0, p- 



_ Po 

3 

P 

d H S 



(24) 



= 0, 



where is the deviator of the Hencky strain tensor. 
The experimental study of a carbon-fibre epoxy compos- 
ite, shocking along the through thickness orientation axis, 
showed no evidence of an inelastic deformation. There- 
fore, elastic constitutive relations (before and after the 
transition zone) are considered here for approximating 
Hugoniot Stress Levels (HSLs) behind the shock wave: 



2p5 



^ijkl £ i'j £ kli a ij —P 



de 

Ik"-' 



1 

2p 



r H r H 

T ij E klJ 



(25) 



r0> ^ijkli a ij> Pij 



A f~iA A nA 

Po j ^ijkh a iji Piji u l 

Pa^ijkli a ijT PijJ u l 



< u. 

> u: 



(26) 

where afj is the Hencky stress tensor corresponding to 
the Hencky strain tensor, p^ is the initial density, C* jfc; is 
the elastic stiffness matrix, a*„- and Bf A are the first and 
second generalisations of the Kronecker delta symbol, Pq 
is the initial density of the anisotropic CFC material, Pq is 
the released density of the isotropic (damaged) CFC ma- 
terial, Cfjki is the anisotropic elastic stiffness matrix cal- 
culated using the properties presented in Tabled ^%jki ^ 



TABLE II. Material properties of damaged Carbon Fibre 
Composite (CFC). 



Parameter Description 



CFC 



po [kg/m 3 ] 
A [GPa] 
p [GPa] 
Cb [m/s] 
C L [m/s] 

a{l, «22, «33 



Density 

Parameter Lame (isotropic case) 
Shear modulus (isotropic case) 
Isotropic bulk sound speed 
Sound speed (isotropic case) 
Tensor a%j (isotropic case) 
Lensor (isotropic case) 



1400.0 

10.434 

0.18 

2745.6 

2777.0 

1.0 

1.0 



the isotropic elastic stiffness matrix calculated using the 



properties presented in Table HH af- and f3^ are the first 
and second generalisations of the Kronecker delta symbol 
for the anisotropic CFC material (see Table HJ), ajj and 
0[a are the first and second generalisations of the Kro- 
necker delta symbol for the isotropic CFC material (see 
Table [TTJ) , and u* is the transition particle velocity, the 
later dependent upon the orientation (hence, u* = u* for 
through the thickness orientation and u* = u*f for along 
the fibre 0° orientation). As a result, the generalised 
Hook law for a two-wave structure has the form 



h _ P n * H 

ij ~ T*^ijfei £ fci' 

Po 



e-H _ ti _ * 

ij ~ ij ij 



H * 
°kl a kl 



(27) 



where is the generalised deviator of the Hencky stress 
tensor. Using experimental data for through thickness 
orientation, the following algorithm is developed for an 
accurate extrapolation of experimental (through thick- 
ness orientation) thermodynamic states, (<r^) p , to high- 
pressure shock Hugoniot states, (Ph) p - 



(H\P (p) n * ( 
K a ij) =— T C y33l 

Po 

p _ Pij (Sid) 



s 



H 



'33J 1 



(pT = - 

(p eos Y 



S. 



H 



P*kl a ll 

(^sDexp ~~ V" 33 

CpT-(pT. 
(4) P 



/ a 33> 



(28) 



(ufY 



33) exp 

EOS\P 



(pT 



Po (u p ) p ' 



where the notation (») p denotes interpolation point (or 
experimental point) p, ( cr 33 SL ) ex p represents the experi- 
mentally measured Hugoniot Stress Levels (HSLs) behind 
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TABLE III. EOS data for CFC used in analysis. 



Parameter Description 



CFC 



4 1 


m/s] 


Velocity curve intercept 


3590.6 


s A 




First slope coefficient 


10.755 


c AA 


[m/s] 


Velocity curve intercept 


3590.6 


s AA 


First slope coefficient 


10.755 


r IA 
c 


[m/s] 


Velocity curve intercept 


2745.7 


C IA 
>->l 




First slope coefficient 


2.9119 


7o 




Griineisen gamma 


0.8500 


(I 




First-order volume correction 


0.5000 


Cl [ 


m/s] 


First anisotropic sound speed 


3599.6 


CI I 


[m/s] 


Second anisotropic sound speed 


2264.2 


[ 


m/s] 


Velocity curve intercept 
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the longitudinal (through thickness orientation) shock 
wave. The data for the bulk shock wave propagation 
Iff (for a single wave and two- wave structures) in the se- 
lected CFC material has been obtained and subsequently 
fitted by straight lines, that is, using linear relationships 
of the form ([29]) and (see Figures Q, ©). 
For a single wave structure: 

U A = c A + Sfu P) c A , S A > 0, V u p . (29) 
For a two-wave structure: 



U 



„AA 



~<AA„ 



-S(-u 
S[ A 



c AA ,s AA >o,u p < u ;, 



V 



JA qIA 



pi H) 



, S( A >0, u p >uZ 



(30) 



The EOS data for the selected CFC material is pre- 
sented in Table IIIII It is important to point out that 
isotropic (damaged composite) has ambient released den- 
sity p 1 = 1.40 g/cm 3 , which is less than the original 
density p A — 1.50 g/cm 3 . This fact is explained by the 
damage softening process. 

Finally, having obtained all of the EOS data in terms of 
the shock and particle velocities, the following algorithm 
is used to obtain an accurate extrapolation of the high- 
pressure shock Hugoniot states to other thermodynamic 



states (HSLs) for the selected shocked carbon-fibre epoxy 
composite (CFC): 



p_ 

Po 



^ijkl e kh 



ft) 



(31) 



U A = c A + S A u p , P A = p*Ufu p , 



and 



EOS 



Cp*) 



-EOS\ 



mn mn 



i + l 



HSL\ 
ij ) 



-l c y 



kl^kl 



(32) 



where (cr^ ) represents the trial Hugoniot Stress Levels 
(HSLs) behind the shock wave, crfj SL represents the true 
Hugoniot Stress Levels (HSLs) behind the shock wave, p 
is the density obtained from (|23p . Pq is the initial den- 
sity (either p A or pg), p is the relative change of volume 
calculated according to (fT5j) and T (v) is the Griineisen 
gamma calculated according to ([15]) . The iterative algo- 
rithm (|32p is performed until the following convergence 
criterion is achieved: (e) f+1 — (e)^ < |(e) i |- error, where 
the notation (•) j and represents, respectively, phys- 

ical quantities in (f3~lj) and (f3"2"| at the i and (i+l) iterative 
steps, meanwhile error = 10 -5 represents the numerical 
error of the iterative algorithm. 

In the following figures, the proposed anisotropic equa- 
tion of state (EOS) for the selected carbon-fibre epoxy 
composite (for a single wave and two-wave structures) is 
examined in terms of the (i) shock and particle velocities 
and (ii) stress (or pressure) and particle velocities. Fig- 
ures ((4]) and ([6]) display, for a single wave and two-wave 
structures, respectively, the relationships between shock 
velocities and particle velocities for (i) through thickness 
and (ii) along the fibre 0° orientation. These figures also 
show the relationship between the anisotropic generalised 
bulk shock wave and the particle velocities. Figures 
([5]) and (O depict, the for a single wave and two- wave 
structures, respectively, the shock Hugoniot states for (i) 
through thickness Pjj (Hugoniot pressure), P£ (total gen- 
eralised pressure), p^os (EOS pressure) (ii) along the 
fibre 0° orientation Hugoniot states (Hugoniot pres- 
sure), Pp (total generalised pressure), Pp OS (EOS pres- 
sure) and, also, the anisotropic generalised bulk shock 
Hugoniot states, Pjj. From Fig. it follows that 

the anisotropic generalised bulk shock wave (for a sin- 
gle wave structure) has a higher velocity of propagation, 
for the selected carbon-fibre epoxy composite, compared 
with the longitudinal (through the thickness) shock wave, 
for all experimental and fitted points. The EOS value of 
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c A in ([29]) was determined to be 3590.6 m/s (for a sin- 
gle wave structure). However, Fig. ([6]) shows different 
CFC behaviour for a two- wave shock structure, where 
the anisotropic generalised bulk shock wave has a lower 
velocity at low particle velocity and a higher velocity at 
high particle velocity. The EOS values (for a two-wave 
structure) of the velocities c AA and c\ A in (|30[) were 
determined to be 3590.6 m/s and 2745.7 m/s, respec- 
tively. In isotropic metals, the empirically derived EOS 
value of Co equates with the theoretical bulk sound speed. 
Note that, for a two- wave structure, the fitted EOS value 
of Cq" 4 = 2745.7 m/s in the isotropic region equates 
with the theoretical bulk sound speed cb — 2745.6 m/s, 
and is lower than the measured anisotropic longitudinal 
(through the thickness) sound speed of 3020 m/s. For the 
anisotropic region, the values of c A and c AA given above 
are significantly greater than the measured longitudinal 
(through the thickness) sound speed of 3020 m/s. These 
values of c A and c AA are also greater than the fitted EOS 
values of Cq = 3228.5 m/s and c F — 3567.7 m/s (for a 
single wave structure) and c AL — 3274.0 m/s (for a two- 
wave structure) - this applies to the longitudinal through 
thickness and along the fibre 0° orientations. However, 
c A and c AA are smaller than the fitted longitudinal (along 
the fibre 0° orientation) EOS value of c AF = 3933.1 m/s 
(for a two- wave structure). Note that the longitudinal 
(through thickness) fitted EOS value of c L = 3230 m/s 
is greater than the measured longitudinal (through thick- 
ness) sound speed of 3020 m/s [23j. This is a behaviour 
that has been observed in many polymers, including 
epoxy resins [HI], [13]. Hence, for a two- wave structure 
containing damage softening effects, similar conclusions 
can be observed in many anisotropic polymers, that is, 
the longitudinal (through thickness) EOS value of c AL 
will be greater than the measured longitudinal (through 
thickness) sound speed. The longitudinal (along the fi- 
bre 0° orientation) fitted EOS values of c F = 3567.7 m/s 
(for a single wave structure) and c AF — 3933.1 m/s (for 
a two- wave structure) are smaller than the respective cal- 
culated longitudinal (along the fibre 0°, 90° orientations) 
sound speeds of 6762.0 m/s and 6666.0 m/s. Further- 
more, the fitted EOS values of c A = c AA = 3590.6 m/s 
may be compared with the first and second generalised 
anisotropic bulk speeds of sound, as the generalisation of 
an isotropic case [25j: 



cq,c aa e [cu,ci], a 



K c K s 

, c/7 = a . (33) 

Po V 



Using the CFC elastic properties presented in Table U it 
can be seen that the fitted EOS value of c A — c AA = 
3590.6 m/s is in the range [c//,cj], where the analytical 
calculations give c/j = 2264.2 m/s and c/ = 3599.6 m/s. 
It is important to re- iterate that, in the isotropic region, 
cjj = cj = 2745.6 m/s meanwhile the fitted EOS value 
of Co~ A = 2745.7 m/s was obtained. 



The experimental shock velocity in the fibre 0° orien- 
tation is initially greater than that corresponding to the 
through thickness orientation. In time, the shock velocity 
decreases with pressure and, eventually, there is conver- 
gence between these data sets. A number of mechanisms 
have been proposed to explain this behaviour [21, 22, 23|. 
However, the experimental data shows that, at lower par- 
ticle velocities, the stress pulse is transmitted through an 
anisotropic mixture of epoxy binder and fibres (see Ta- 
ble [l|, whereas at higher particle velocities, this pulse 
is transmitted through an isotropic mixture instead (see 
Table [H} . The stable shock waves (where shock velocity 
decreases with pressure) can exist when the shock front 
breaks up into two or more waves 30] . As is shown above, 
the two- wave front structure is sufficient to fit experimen- 
tal data for the orientations through the thickness and 
along the fibre 0°. 

Using a generalised decomposition of the stress ten- 
sor, the generalised anisotropic bulk Hugoniot pressure 
(Pjj = PqII a u p ) can be defined and compared to the lon- 
gitudinal Hugoniot pressures (Pjj = p^UgUp for through 
the thickness; Pjj = PqII f u p for along the fibre 0° orien- 
tation), as shown in Fig. ([5]) for a single wave structure 
and in Fig. jJJ for a two-wave structure. It can be seen 
that, for a single wave structure, there is a significant 
difference between the longitudinal Hugoniot pressures, 
Pjj and Pjj, and the calculated generalised anisotropic 
bulk Hugoniot pressure, Pjj. This indicates that, for 
highly anisotropic materials and the assumption of a 
single wave structure, the anisotropic bulk shock front 
will be supersonic with respect to the longitudinal shock 
front in the least stiff direction (e.g., the through thick- 
ness orientation). However, this is not the case in the 
two-wave structure. According to the generalised de- 
composition (TT]), the stress tensor has been split into the 
generalised spherical component, p*ctij, and the gener- 
alised deviatoric component, Sij. The total generalised 
anisotropic "pressure", p* , comprises a sum of two terms, 
p* = p EOS + p s , where p EOS corresponds to the thermo- 
dynamic (EOS) response calculated from (fT4"]) - (fT5"|) and 
p corresponds to the generalised deviatoric stress calcu- 
lated using (J5]). 

Figures (JSJ and © display p* and p EOS . It can be 
seen that there is an increasing divergence of p* from 
pEOS j ue to g rea t er p ar ticle velocities associated with 
an increased contribution from p . In addition, it can be 
seen that pressure p EOS for the selected CFC material 
is greater than the generalised total pressure, p* . Fur- 
ther comparison of different pressures (see Figs. ([5]) and 
([7]) ) shows that there is an increasing divergence between 
the measured Hugoniot stress (HSLs) and the calculated 
pressures Pjj, Pjj, p FOS ', pf os ', p* F and p* L - this cor- 
responds to the increasing contribution from the gener- 
alised deviator of the stress tensor, Sy . In addition, it is 
important to analyse the difference between the shock ve- 
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FIG. 4. Anisotropic generalised bulk U a shock velocity, lon- 
gitudinal (through the thickness) Ug shock velocity and lon- 
gitudinal (fibre 0° orientation) U f shock velocity for a single 
wave structure of carbon fibre-epoxy composite under shock 



loading, where U s is defined by (|29 



locitics in the two-wave structure at the phase transition 
point, these evaluated using: 
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FIG. 5. Comparison of the different Shock Hugoniot pres- 
sures Pk, Pt, Pf OS , Pf , Pf, Pf° S and P£ for a single- 
wave structure of the carbon-fibre epoxy composite behind the 
shock wave in stress-particle velocity space. The Pjj, P^ and 
PjJ curves are calculated from the equations P^ = poUgUp, 
Px = PqUs u p and P^ = PqU^u p respectively, where 
is defined by fl29j. The P£, Pf os , P F and P# os curves are 
defined by JTJJ) and (Ti3)) . 



[ttF~I ( AF IF\ | I Q AF qIF\ 
[U s \ = (Co - c J + [ S l ~ S l ) ' 



(35) 



where [C//'] tl _ u * = 154.935 rajs for the through thick- 
ness orientation and \\Jz\ , = 754.705 m/s for the 

fibre 0° orientation at the transition points and u*^ 
respectively. Using fitted EOS data (see Table [TTT|) . the 
difference between shock waves at the transition point 
for the through thickness orientation is smaller than that 
along the fibre 0° orientation. Therefore, the material 
behaviour through the thickness orientation can be ap- 
proximated with a single wave structure (see Figure [4]), 
this confirmed during the examination of Hugoniot Stress 
Levels (HSLs) through the_thickness orientation in terms 
of the particle velocities [33]. 

Finally, in Figs. (JSJ) and ©, the effect of orientation 
on the Hugoniot Stress Levels (HSLs) in stress-particle 
velocity space is examined for single wave and two-wave 
structures. It can be seen from the experimental points in 
Figs. (j5J), © that, at lower stresses, the mixture of epoxy 
binder and 0° fibres is stiffer, however, as stress increases, 
the experimental Hugoniots for a mixture (CFC compos- 
ite) of both orientations converge. This is in agreement 
with the behaviour of the shock velocities shown in Fig. 
© for a two-wave structure. It is also clear from Fig. 
© that a single wave structure is not capable of predict- 
ing correctly the Hugoniot Stress Levels (HSLs) along 



the fibre 0° orientation to agree with the stability re- 
quirements formulated by Bethe (30j . It is important to 
note that there is no true prediction of the experimental 
data shown in figures and (J9j) for through the thick- 
ness orientation. These experimental points (through the 
thickness) have been used to define the material param- 
eters in the presented anisotropic EOS model. There is 
only a true prediction of the experimental data shown 
in figures © and for along the fibre 0° orientation. 

Figures ©, ^ show qualitatively that the anisotropy 
of a composite material (carbon fibre-epoxy composite) 
has a strong effect on the accurate extrapolation of high- 
pressure shock Hugoniot states to other thermodynamic 
states for shocked anisotropic composite (CFC) materials 
(e.g., a carbon-fibre epoxy composite) of any symmetry. 



CONCLUSIONS 

An anisotropic equation of state is proposed for the 
accurate extrapolation of high-pressure shock Hugoniot 
states to other thermodynamic states, for a shocked 
carbon-fibre epoxy composite (CFC) of any symmetry. 
The proposed equation of state, which uses a generalised 
decomposition of the stress tensor [3, H^, \ 



27, 2t 



represents the mathematical and physical generalisation 
of the Mie-Griineisen equation of state for an isotropic 
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FIG. 6. Anisotropic generalised bulk U a shock velocity, lon- 
gitudinal (through the thickness) U^ shock velocity and lon- 
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gitudinal (fibre 0° orientation) U s shock velocity for a two- FIG. 7. Comparison of the different Shock Hugoniot pres- 



wave structure of carbon fibre-epoxy composite under shock 
loading, where Ug is defined by (f30j 



material, and reduces to the latter in the limit of isotropy. 

Further insights into the anisotropic CFC response un- 
der shock loading can be gained from an examination of 
the material EOS in terms of: shock stresses (total gen- 
eralised anisotropic pressure, p*; generalised anisotropic 
bulk Hugoniot pressure, P^; pressure, p EOS , correspond- 
ing to the thermodynamic (equation of state) response; 
pressure, p s , corresponding to the generalised deviatoric 
stress), shock velocities (shock velocity in the through- 
thickness orientation, Ug, shock velocity in the fibre 0° 
orientation, Ug , and the generalised anisotropic bulk 
shock velocity, U A ) and the particle velocity, u p . Fig- 
ure ([6]) shows linear relationships, for a two-wave struc- 
ture, between the shock velocities ,Ug ,U A and the 
particle velocities, u p , over the range of measurements 
made during experiments. The values of c AL ,c AF ,c AA 
and Cq L ,Cq F ,Cq A were determined. The values c AA and 
Cq A (the intercept of the Uf-u p curve for two-wave struc- 
ture) are in the interval between the first and second 
generalised anisotropic bulk speeds of sound [25| (for 
non-linear anisotropic elastic and isotropic elastic shock 
waves). This is a behaviour that is observed in many 
polymers, including epoxy resins. When c AL , c AF and 
c AA are compared to the measured longitudinal sound 
speed Cl (Table HI Table |H*I)) . it can be seen that the 
former values are significantly greater than the latter, in 
the through thickness orientation. This indicates that, 
for highly anisotropic materials, anisotropic shock fronts 
(at lower particle velocity) are always supersonic with 
respect to the longitudinal sound speed in the least stiff 
direction (for example, longitudinal sound speed in the 
through thickness orientation). It is possible that the 



sures P%, P£, Pf os , Pg, P F , P# os and P§ for a two- 
wave structure of the carbon-fibre epoxy composite behind 
the shock wave in stress-particle velocity space. The Pf}, P^ 
and P£ curves are calculated from equations Pj} = p%U^u v , 

p respectively, where U A is 
, Pi and Pp OS curves are 



P'h = plU'up and P£ = p* U, 



defined by (001 ■ The P£, 
defined by (T2} and 
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FIG. 8. Hugoniot Stress Levels (HSLs) for a carbon-fibre 
epoxy composite behind the shock wave in the through thick- 
ness and fibre 0° orientations in stress-particle velocity space, 
assuming a single wave structure. The experimental data for 
Hugoniot Stress Levels (HSLs) was obtained by Millett et al. 
[23J . The dotted curve is calculated using the anisotropic EOS 
proposed in (|31|l , (|32[) and experimental data for the through 
thickness orientation Ug — u p . The solid curve is calculated 
using the anisotropic EOS proposed in (|31|l . (132p and experi- 
mental data for the fibre 0° orientation Ua — u v . 
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FIG. 9. Hugoniot Stress Levels (HSLs) for a carbon-fibre 
epoxy composite behind the shock wave in the through thick- 
ness and fibre 0° orientations in stress-particle velocity space, 
assuming a two-wave structure. The experimental data for 
Hugoniot Stress Levels (HSLs) was obtained by Millett et al. 
[23| . The dotted curve is calculated using the anisotropic EOS 
proposed in (|31[l , (|32[) and experimental data for the through 
thickness orientation U§ — u p . The solid curve is calculated 
using the anisotropic EOS proposed in (f3Tj) . (f32|) and experi- 
mental data for the fibre 0° orientation U f — u p . 

generalised anisotropic bulk shock velocity, Uf, depends 
non-linearly on particle velocity for a given anisotropic 
material (or composite materials) , as has been shown for 
some isotropic polymers such as PMMA [l[ - unfortu- 
nately, the corresponding experimental data for carbon 
fibre materials (or other anisotropic materials) cannot be 
located at this time. 

An analytical calculation showed that the Hugoniot 
Stress Levels (HSLs) for a carbon-fibre epoxy compos- 
ite do not agree with the experimental data for a sin- 
gle wave structure methodology for the stability require- 
ments formulated by Bethe [30j] . However, the material 
behaviour in the through-thickness orientation can be ap- 
proximated by a single wave structure due to the lack of 
any significant discontinuity at the transition zone [33j |. 
This approximation will not take into account changes in 
the material elastic properties during the damage soft- 
ening process. In addition, an analytical calculation 
showed that the Hugoniot Stress Levels (HSLs) in dif- 
ferent directions, for a CFC composite subject to the 
two-wave structure (non-linear anisotropic and isotropic 
elastic waves) , agree with experimental measurements at 
both low shock intensities (where the 0° orientation was 
significantly stiffer than the through-thickness orienta- 
tion) and at high shock intensities (where the HSLs of 
the two orientations converged due to the presence of 
damage softening), this also in agreement with the sta- 



bility requirements formulated by Bethe [3C 
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